Theory of plasmon-enhanced high-harmonic generation in the vicinity of metal 

nanostructures in noble gases 
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We present a semiclassical model for plasmon-enhanced high-harmonic generation (HHG) in the 
vicinity of metal nanostructures. We show that both the inhomogeneity of the enhanced local fields 
and electron absorption by the metal surface play an important role in the HHG process and lead 
to the generation of even harmonics and to a significantly increased cutoff. For the examples of 
silver-coated nanocones and bowtie antennas we predict that the required intensity reduces by up 
to three orders of magnitudes and the HHG cutoff increases by more than a factor of two. The 
study of the enhanced high-harmonic generation is connected with a finite-element simulation of the 
electric field enhancement due to the excitation of the plasmonic modes. 

PACS numbers: 42.65.Ky, 78.67.Bf 



Progress of contemporary ultrafast laser physics has 
led to the generation of attosecond laser pulses by high- 
harmonic generation (HHG) [lj opening the way to ex- 
citing opportunities to study time-resolved fundamental 
electronic processes in atoms and molecules, and to in- 
vestigate the collective and correlated dynamics of many- 
electron systems on the attosecond scale On the 
other hand, in an another rapidly evolving field metallic 
nanostructures can be used to achieve localization of light 
on sub-wavelength nanometer scale and to realize huge 
plasmonic field enhancements in the vicinity of these ob- 
jects. At present plasmonics has been the subject of ex- 
tensive theoretical and experimental research and experi- 
ences a dramatic growth of applications 0-Gjl. Recently 
an interesting experiment has been reported based on a 
connection of both fields: high-order harmonic genera- 
tion (HHG) has been demonstrated by nJ pulses directly 
from a laser oscillator by exploiting the local field en- 
hancement near a metallic bowtie-shaped gold nanoele- 
ment [13l ] . The elimination of complex and expensive 
amplifiers and the opportunity to perform HHG and at- 
tosecond experiments at MHz repetition rates can sig- 
nificantly extend the ability of extreme VUV genera- 
tion and attosecond pump-probe spectroscopy and will 
lead to novel applications in lithography and imaging. 
However, a theoretical description of plasmon-enhanced 
HHG, which is a key for further progress in this field, is 
currently missing. 

In this paper we present a theoretical investigation of 
plasmon-enhanced HHG in a noble gas in the vicinity 
of a silver nanocone and a silver bowtie structure. Our 
goal is directed to an understanding of the role of metal 
nanostructures and the inhomogeneous plasmonic field 
enhancement in the HHG process. We describe the mi- 
croscopic laser-atom interaction by the time-dependent 
dipole moment in the semi-classical approximation 14 1. 



In the vicinity of the metal nanostructures the input laser 
intensity is enhanced by up to three orders of magnitude 
by the plasmons in so-called "hot spots" and sensitively 



depends on the position. This strongly enhanced field 
ionizes the noble gas atoms, the resulting electrons are 
accelerated in the enhanced field and recollide with the 
parent ion, transforming their kinetic energy into high 
harmonic radiation. In the hot spot, electron experience 
the inhomogeneous electric field, and some of them will 
hit the metal surface and never return to the parent ion. 
We show that both field inhomogeneity and the electron 
absorption at the metal surface lead to the emission of 
even harmonics as well as to more than a twofold increase 
of the cut-off frequency. 

We describe HHG in the framework of a model that 
is an extension of the Lewenstein model modified to ac- 
count for the field inhomogeneity in the vicinity of the 
nanostructure, as well as for the electron absorption by 
the metallic surface. Under the action of the electric 
field E(t), the electron leaves the atom at the time mo- 
ment t s and later recombines with the parent ion at the 
time tf. The maximum HHG contribution is provided 
by electrons with the canonical momentum p for which 
their position at tf coincides with that of the parent ion. 
Neglecting the Coulomb potential and using the station- 
ary phase method the time-dependent high- harmonic 
dipole moment is given by 
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Here e is an arbitrary small parameter, w is the cen- 
tral frequency, d x (p) = i2 7 - 25 [huJomll p ] 5 / 4: ir~ 1 p/(p 2 + a) 3 

with a = 2ml^ 2 Ip^ 2 is the dimensionless dipole matrix 
element, and S(t f ,t s ) = f*J(I p + (p(t) - eA{t)) 2 /2m e )dt 
is the classical action equal in the Lewenstein model to 
S (tf,t s ) = I p At - 0.5e 2 (AB 2 /Ai + AC)/m e . Here 
A = dA(t)/dt = E(t), B(t) = A(t), C(t) = A 2 (t), for any 
function F we define AF = F(t f ) - F(t s ), At = tf-t s , 



2 



and I p is the ionization potential. For the unmodi- 
fied Lewenstein model, in the stationary-phase approx- 
imation the canonical momentum is given by eAB / At. 
H(tf,t s ) is a function describing the electron absorption 
at the metal surface; in the unmodified Lewenstein model 
we have H(tf,t s ) = 1. We extend the semiclassical HHG 
model accounting for the inhomogeneity of the field in the 
hot spot by adding the first-derivative term to the field 
E(t,x) = E(t)(l + x/di n h), and consider this derivative 
term as a perturbation, including only the first-order cor- 
rection to the electron trajectory. As a result, the expres- 
sions for the momentum p(t s ) and S(tt,t s ) are modified 
as follows: 

p(t s ) = e[A(t s ) + 

AB - A(t s )At + /3{0.5(AB) 2 - AD + C(t s )At) 

At - P{2AG + At[B{t) + B{t s )]) J ' ( ' 

S(t f ,t s ) = S (t fl t s ) + e 2 /3m; 1 {p{t s ) 2 [2AG 

+ (B{t f ) + B{t s ))At] + p{t s )[-ACAt 
+2AD] + (C(t f ) + C{t s ))AB - 2AF}, (3) 

where D(t) = C(t), F{t) = C{t)E{t), G{t) = B{t), 
j3 = e/ (m e di n h), and To is the optical period. Another 
modification of the model is described by the function 
H(t,t s ) under the integral in Eq. (1). This function is 
equal to 1 unless the electron during the motion hits 
the metal surface positioned at d sur , otherwise we as- 
sume that the electron is absorbed by the surface and set 
H(t,t s ) = 0. We have used the finite-element Maxwell 
solver JCMwave to model the plasmonic field enhance- 
ment in the vicinity of the nanostructure, taking the ex- 
perimental complex-valued dielectric function of silver. 
For short input pulses, to include the dispersion of silver 
and of the field distribution near the nanostructure, we 
calculated the field distribution for a set of wavelengths 
within the input pulse spectrum and then reconstructed 
the position-dependent pulse profile from the wavelength- 
and position-dependent field. In this paper, we disregard 
the modification of the field induced by the neighboring 
structure, the depolarization effects are weak in our case. 
To prove that the mesh size is small enough, we halved 
it and found no significant change in the results. 

First we analyze the principal effects described by the 
above extended model in dependence on the model pa- 
rameters di n h and d sur . To understand the effect of the 
inhomogeneity of the field and of the metal surface on 
the HHG process, in Fig. 1 we consider the HHG spec- 
tra for different parameters di n u [Fig. 1(a)] and d sur [Fig. 
1(b)] for cw radiation with a typical intensity in the hot 
spot of 200 TW/cm 2 . In Fig. 1(a) even for a very low 
inhomogeneity, characterized by di n h — 1000 nm, the in- 
verse spatial symmetry is broken, which leads to a qual- 
itatively new effect: the generation of even harmonics 
(shown by the blue crosses). For larger inhomogeneities 



with di n h = 20 nm the amplitude of the even harmon- 
ics approaches that of odd harmonics. Simultaneously, 
the harmonic cutoff shifts to higher harmonic numbers 
and become much less pronounced. Due to inhomogene- 
ity, electrons have higher velocities when they reach the 
parent ion, leading to a higher kinetic energy and there- 
fore to a higher cut-off frequency. The influence of the 
metal surface and the associated electron absorption is 
studied in Fig. 1(b). One can see that for a distance to 
the surface d sur — 3 nm, no electron reaches the surface, 
and correspondingly no modification in the dipole mo- 
ment and no even harmonics appear. For shorter d sur , 
some electrons will reach the surface. Since the function 
H(t,t s ) introduces a sharp modulation in the temporal 
profile of the harmonics, an extension of the spectrum to 
higher harmonic numbers can be seen in Fig. 1(b), com- 
bined with the emission of even harmonics. From Fig. 1 
one can see that the combined effect of the field inhomo- 
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FIG. 1: Influence of the field inhomogeneity (a) and of the 
metal surface (b) on the HHG process. Odd harmonics 
are presented by curves while even harmonics are shown by 
points. In (a), the field has the inhomogeneity scale d in h of 
1000 nm (red solid curve and red vertical crosses), 20 nm 
(green short-dashed curve and green diagonal crosses), and 
5 nm (blue long-dashed curve and blue squares). In (b), 
the distance to the metallic surface d sur is 3 nm (red solid 
curve), 2.62 nm (green short-dashed curve and green diag- 
onal crosses), and 1.5 nm (blue long-dashed curve and blue 
squares). For both cases, cw radiation at 830 nm with inten- 
sity in the hot spot of 200 TW/cm 2 in argon is considered. 
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FIG. 2: High harmonic generation in the vicinity of a 
nanocone structure. In (a), the scheme of the nanocone is 
shown. The curvature radii and the thickness of silver coat- 
ing are 5 nm, the y-polarized field propagates along the x 
direction. In (b) the intensity enhancement is presented for 
the zy cross-section of the incident field at 830 nm. In (c) 
and (d), the distribution of the harmonic cutoff is shown for 
the model without (c) and with (d) modifications by the field 
inhomogeneity and metal surface for 10-fs input pulses with 
0.3 TW/cm 2 intensity and argon surrounding the tip. In (e) 
the cutoff for the modified (green curve) and unmodified (red 
curve) model is shown along the y axis, and in (e) the spec- 
trum for cw excitation with the above parameters is shown. 



geneity and of the metal surface can increase the HHG 
cutoff by a factor larger than two. 

Next we perform a realistic simulation of HHG in the 
vicinity of a silver nanocone as depicted in Fig. 2(a). 
Such nanocone structure can be produced by silver de- 
position on a chemically created silicon needle without 
the need of nanoscale lithographic control. The curva- 



ture radius of the nanostructure edges is the key param- 
eter determining the maximum field enhancement, with 
optimum structure having the smallest curvature radius. 
We have assumed radii of 5 nm, in correspondence with 
the current state of the manufacturing technique. The 
spatial field distribution shown in Fig. 2(b) with pa- 
rameters given in the caption demonstrates intensity en- 
hancement factors of up to 1000 in the "hot spot" in 
the vicinity of the tip. From the comparison of xy and 
zy cross-sections (not shown) one can see that the field 
enhancement is higher on the surface directed towards 
the incoming beam, but this is a minor effect compared 
to the strong enhancement near the tip. Note that the 
field is inhomogeneous on the scale below 10 nm in the 
"hot spot" which is located directly above the metal sur- 
face. The parameters di n h{x,y, z) and d sur {x,y, z) are 
calculated from this spatial distribution of the field and 
are used for the HHG simulation by using the extended 
model as given by Eq. (l)-(3). In Fig. 2(c) and (d) 
we present the spatial distribution of the harmonic cut- 
off N t h for an excitation by 830-nm sech-shaped pump 
pulses with intensity of only 0.3 TW/cm 2 , FWHM of 
10 fs, and polarization and wavevector as indicated in 
Fig. 2(a). The field enhancement in the "hot spot" 
increases the intensity to roughly 300 TW/cm 2 , which 
for the chosen wavelength corresponds in the unmodified 
Lewenstein-model to a cut-off of roughly the 50 th har- 
monic, as can be seen in Fig. 2(c). However, due to 
the effect of the field inhomogeneity and the presence of 
the metal surface, the cutoff increases up to 105, which 
correspond to the emission of 150-eV photons using an 
input intensity of only 0.3 TW/cm 2 . Such pulses can 
be obtained from a typical laser oscillator without am- 
plifier. In Fig. 2(e) a comparison between the extended 
model including the effects of the field inhomogeneity and 
metal surface, and the unmodified Lewenstein model is 
drawn. The cut-off harmonic number as a function of 
the position reaches values of above 100 for the mod- 
ified model, while the unmodified only predicts values 
around 40. Further from the tip of the nanocone, the 
field become relatively homogeneous, and the distance to 
the metal surface is large, therefore both models coin- 
cide. Note however, that although the modifications of 
the Lewenstein model are important only in a small frac- 
tion of the volume, it is exactly this volume ("hot spot") 
which dominate the generation of high harmonics and 
therefore these modifications become especially impor- 
tant. This is clearly visible in Fig. 2(f), where we show 
the spectrum of the high harmonics emitted in the x di- 
rection integrated over the nanostructure position, where 
a significant emission of even harmonics as well as the ex- 
tended and smoothed cutoff can be seen. We have also 
calculated the phases of the HHG harmonics separately 
for short and long elelctron trajectories, and found that 
the modification of phases due to the field inhomogeneity 
is negligible for short trajectories but reaches 7r/2 for the 
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FIG. 3: Temporal shape (a) and the temporal shape after 
filtering out low harmonics (N <40) of the output(b) with 
input intensity of 0.3 TW/cm 2 and duration of 10-fs at 830 
nm for a nanocone structure. 



long ones. This, however, will not influence the diver- 
gence of the HHG radiation, since the hot spot is smaller 
than the typically considered wavelength and is a point- 



like source. In [15( we calculated the variation of the 
high-harmonic phases with different intensities depend- 
ing on the transverse positions in the beam. It is shown 
that in the considered example the phases change only 
slightly for different intensities within the beam for har- 
monic numbers below 50, which can be compensated by 
an appropriate curvature of the pump beam. For higher 
harmonic numbers the phases show strong changes, re- 
sulting in beam distortion due to irregular spatial mod- 
ulation. 

In Fig. 3, the temporal profile of the output harmonic 
radiation is shown for 10-fs input pulses. The tempo- 
ral profile in Fig. 3(a) is dominated by low-order (<10) 
harmonics. It is well-known that filtering out the lower 
harmonics leaves only few electron trajectories with max- 
imum kinetic energy at the recollision time, which can 
lead to the formation of attosecond pulses. In Fig. 3(b) 
we present the temporal profile for the above case but 
with harmonics below the 50th order filtered out. One 
can see the formation of few pulses with down to 150- 
attosecond FWHM. These pulses are shorter and more 
irregular than those obtained from the unmodified model 
(not shown). 

In order to compare the results of our model with ex- 
perimental measurements, we have simulated high har- 
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FIG. 4: High harmonic generation in the vicinity of a bowtie 
structure. In (a), the scheme of the bowtie structure and its 
geometric parameters are shown, the :r-polarized field propa- 
gates along y direction. In (b) the field enhancement is pre- 
sented for the xy cross-section of the incident field at 830 nm; 
in (c), the distribution of the harmonic cutoff is shown for 10- 
fs input pulses with 0.5 TW/cm 2 intensity for HHG in argon 
surrounding the bowtie structure. 



monic generation in the vicinity of a bowtie nanoantenna, 
such as used in Ref . fl~3j . Structure and parameters of the 
studied bowtie antennas is shown in Fig. 4(a) consisting 
of two silver isosceles triangles with curvature radius and 
distance between the vortices in the range of 10 nm. In 
Fig. 4(b), the field enhancement by the bowtie antenna 
with parameters as given in the figure is presented, which 
shows a " hot spot" between the inner vortices of the tri- 
angles and a maximum enhancement of about 10 2 . For 
other geometrical parameters of the structure, the en- 
hancements drops with increasing curvature radius and 
distance between the triangle tips. By using this field 
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enhancement HHG is calculated by using the extended 
model. In Fig. 4(c) for incident 10-fs pulses with inten- 
sity of 0.5 TW/cm 2 the spatial distribution of the thresh- 
old harmonic number N t ^ is presented. The calculated 
maximum harmonic number is 23, which corresponds to 
the wavelength of 36 nm. This is in good agreement 
with the experimental findings of [l3j |. where harmonics 
down to 45 nm were observed. Note that both the ge- 
ometry of the bowtie antenna and the input intensity are 
not exactly defined in [l3| which can lead to deviations 
between theory and the experiment. In contrast to the 
nanocone, in the case of the bowtie antenna a relatively 
low field enhancement due to larger curvature radii leads 
to a smaller amplitude of electron motion which, together 
with a relatively homogeneous field distribution, makes 
the model modifications less important than in the case of 
the considered nanocone. Our calculations predict that 
the cutoff value has increased by only roughly 20% due to 
the effects of the inhomogeneity and the metallic surface. 



[1] M. Hentschel, R. Kienberger, C. Spielmann, G. A. Rei- 
der, N. Milosevic, T. Brabec, P. Corkum, U. Heinzmann, 
M. Drescher, and F. Krausz, Nature 414, 509 (2001). 

[2] A. Paul et al, Nature 421, 51 (2003). 

[3] X. Zhang et al, Nat. Phys. 3, 230 (2007). 

[4] P. B. Corkum and F. Krausz, Nat. Phys. 3, 381 (2007). 

[5] H. Kapteyn, O. Cohen, I. Christov, and M. Murnane, 
Science 317, 775 (2007). 

[6] M. F. Kling and M. J. J. Vrakking, Annu. Rev. Phys. 
Chem. 59, 463 (2008). 

[7] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 
(2009). 

[8] G. Sansone et al, Nature 465, 763 (2010). 

[9] S. A. Maier, "Plasmonics: Fundamentals and applica- 
tions" , Springer , Berlin 2007. 
[10] M. Pelton, J. Aizpurua, G. Bryant, Laser Photon. Rev. 
2, 136 (2008). 

[11] S. Nie and S. R. Emory, Science 275, 1102 (1997). 
[12] E. Prodan et al, Science 302, 419 (2003). 
[13] S. Kim et al, Nature 453, 757 (2008). 
[14] M. Lewenstein et al, Phys. Rev. A 49, 2117 (1994). 
[15] Supplementary material is available for this paper at 
|http://prl.aps.org/ 



* Electronic address: gusakov@mbi-berlin.de 



